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Abstract 

At the solar equinox in August 2009, the Composite Infrared Spectrometer 
(CIRS) onboard Cassini showed the lowest Saturn’s ring temperatures ever 
observed. Detailed radiative transfer models show that the observed equinox 
temperatures of Saturn’s A ring are much higher than model predictions as 
long as only the flux from Saturn is taken into account. In addition, the 
post-equinox temperatures are lower than the pre-equinox temperatures at 
the same absolute solar elevation angle. These facts indicate that the A 
ring was not completely cooled down at the equinox and that it is possible 
to give constraints on the size and seasonal thermal inertia of ring particles 
using seasonal temperature variations around the equinox. We develop a 
simple seasonal model for ring temperatures and hrst assume that the internal 
density and the thermal inertia of a ring particle are uniform with depth. The 
particle size is estimated to be 1-2 m. The seasonal thermal inertia is found to 
be 30-50 Jm“^K“^s“^/^ in the middle A ring whereas it is ~ 10 Jm“^K“^s“^/^ 
or as low as the diurnal thermal inertia in the inner and outermost regions of 
the A ring. An additional internal structure model, in which a particle has a 
high density core surrounded by a fluffy regolith mantle, shows that the core 
radius relative to the particle radius is about 0.9 for the middle A ring and 
is much less for the inner and outer regions of the A ring. This means that 


* Corresponding author 

Email address: Ryuji.Morishima@jpl.nasa.gov (Ryuji Morishima) 


Preprint submitted to Elsevier 


September 1, 2015 




the radial variation of the internal density of ring particles exists across the 
A ring. Some mechanisms may be conhning dense particles in the middle 
A ring against viscous diffusion. Alternatively, the (middle) A ring might 
have recently formed (< 10® yr) by destruction of an icy satellite, so that 
dense particles have not yet diffused over the A ring and regolith mantles 
of particles have not grown thick. Our model results also indicate that the 
composition of the core is predominantly water ice, not rock. 
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1. Introduction 

Saturn’s main rings consist of many small particles orbiting around the 
planet. The range of particle sizes deduced from radio and stellar occultations 
is roughly 1 cm to 10 m (Zebker et ah, 1985; French and Nicholson, 2000; 
Cuzzi et ah, 2009). The composition of ring particles is mostly crystalline 
water ice and the mass fraction of contaminants (e.g., Tholins, PAHs, or 
nanohematite) is 10% at most (Epstein et ah, 1984; Cuzzi et ah, 2009) and 
probably less than 1 % (Poulet et ah, 2003). A favorable origin of Saturn’s 
rings with such a high content of water ice is stripping of the icy mantle of a 
Titan-sized satellite (Canup, 2010). This hypothesis indicates that Saturn’s 
rings are primordial (as old as Saturn itself) and the initial mass was a few 
order of magnitude more massive than the present ring mass. As a byproduct 
of viscous diffusion of Saturn’s rings, accretion of icy satellites occurs at the 
outer edge of the rings, followed by orbital expansions of the satellites due to 
the torques induced by the rings and the planet (Charnoz et ah, 2010, 2011). 

While the ring system as a whole may be primordial, local structures of 
the rings generally change on much shorter timescales. Particularly, several 
observational facts indicates that the outermost main ring, the A ring, is 
dynamically young (see Charnoz et ah, 2009 for the review). The A ring 
extends from 122,000 km to 137,000 km from the center of Saturn. Gravita¬ 
tional torques exerted at resonances in the A ring with inner icy satellites, 
such as Mimas, Janus, and Epimetheus, cause the A ring to collapse to its 
inner edge in ~ 10® years (Goldreich and Tremain, 1982; Dones, 1991), unless 
angular momentum is supplied from some other sources. The A ring viscosity 
measured from density wave patterns at satellite resonances is of the order 
of 100 cm^ s“^ or higher (Esposito et ah, 1983; Tiscareno et ah, 2007). This 
also gives a radial diffusion timescale of ~ 10® years for the A ring width (see 
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also Salmon et al., 2010). The high viscosity of the A ring is likely a result of 
wakes produced by gravitational agglomeration of ring particles. Self-gravity 
wakes are observationally suggested from azimuthal brightness asymmetry 
(e.g., Dones et ah, 1993; French et ah, 2007; Colwell et ah, 2006; Hedman et 
ah, 2007; Ferrari et ah, 2009) and also predicted from A^-body simulations 
(Salo, 1995; Daisaka et ah, 2001; Salo et ah, 2004; Robins et ah, 2010). The 
apparently young A ring may or may not be compatible with the primordial 
origin of the ring system. To clarify the origin and evolution of the ring 
system, additional observational constraints are invaluable. 

Since Saturn Orbit Insertion of Cassini in 2004 to the present, the Cassini 
spacecraft has been observing Saturn’s rings in broad wavelengths from ul¬ 
traviolet to radio. The Cassini Composite infrared Spectrometer (CIRS) has 
acquired millions of ring spectra in mid- to far-infrared wavelengths (7 /rm - 1 
mm) at various observational geometries (e.g., Spilker et ah, 2006; Altobelli 
et ah, 2008, 2014; Leyrat et ah, 2008b). The wavelength range of CIRS cov¬ 
ers the Planck peak of the thermal emission, and the effective temperature 
and the geometric Filing factor of a ring are derived by applying a Planck £t 
to a spectrum (see Spilker et ah, 2006 for details). The typical error in the 
derived temperature is much less than 1 K. 

Ring temperatures depend on observational geometric parameters. Ring 
temperatures vary seasonally primarily due to the change of solar elevation 
angle from -27° to 27° and secondarily due to the change of heliocentric 
distance of Saturn (Froidevaux, 1981; Flandes et ah, 2010; Pilorz et ah, 
2014). The face illuminated by the Sun (the lit face) is warmer than the unlit 
face and the difference between the lit and unlit face temperatures increases 
with ring optical depth (Spilker et ah, 2006). Ring temperatures also vary 
diurnally due to eclipse cooling in Saturn’s shadow (Leyrat et ah, 2008b). 
On top of the eclipse cooling, the temperatures of the A ring show quadruple 
azimuthal modulations due to self-gravity wakes (Leyrat et ah, 2008b). Ring 
temperatures increase with decreasing solar phase angle (Altobelli et ah, 
2007; 2009), indicating that ring particles are not isotropic emitters. 

With each temperature variation, different types of constraints on ring 
structures and particle properties are obtained. For example, the diurnal 
thermal inertia of ring particles is estimated to be 10 Jm ^s for all 
main rings, using diurnal temperature variation curves including tempera¬ 
tures in Saturn’s shadow (Ferrari et ah, 2005; Leyrat et ah, 2008a; Morishima 
et ah, 2011, 2014). This low thermal inertia value indicates that surface re- 
golith of ring particles is very fluffy at least for the thickness comparable 
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to the diurnal thermal skin depth, ~ 1 mm. However, fluffiness of deep 
interiors of ring particles cannot be constrained from diurnal temperature 
variations. Very porous regolith on the particle surface is also indicated from 
the opposition effect seen in photometric phase curves (Dean, 2015). 

CIRS radial scans measured at various observational geometries give a 
radial profile of the bolometric Bond albedo (Morishima et ah, 2010). 
The value of is found to be correlated with the ring optical depth, r: 
Ab = 0.1 — 0.4 for the C ring (r ~ 0.1), Ab = 0.5 for the A ring (r ~ 0.5), 
and Ab = 0.6 — 0.7 for the B ring (r > 0.7). This correlation probably 
snggests that Saturn’s rings have been continnously polluted by meteoroid 
bombardments and that pollntion is more effective for optically thinner rings 
(Cuzzi and Estrada, 1998; Elliot and Esposito, 2011). The radial variation 
of Ab inside the A ring is fonnd to be qnite small (Morishima et ah, 2010). 
This is consistent with the small radial variation of visual and near-infrared 
reflectances of the A ring (Porco et ah, 2005; Hedman et ah, 2013). 

At the solar eqninox in Angust 2009, Satnrn’s rings revealed the lowest 
temperatnres ever observed; there was no temperature difference between 
the sonthern and northern faces (Spilker et ah, 2013). At the eqninox, as 
the Satnrn flnx is dominant, the ring temperatnre decreases with increasing 
satnrnocentric distance. The eqninox temperatnre is fonnd to be relatively 
higher at optically thinner regions, because ring particles in thin rings can 
be heated by both the sonthern and northern hemispheres of Satnrn. The 
observed eqninox temperatnres can be well reprodnced by the models used 
in Spilker et ah (2013), except for the A ring. They applied a mnlti-particle- 
layer model developed by Morishima et ah (2009) to the eqninox data of the 
A ring and fonnd that the observed A ring temperatnres are mnch higher 
than model predictions in which only the energy sonrce at the eqninox is as¬ 
sumed to be the thermal and solar-reflected fluxes from Saturn. The eqninox 
temperature anomaly is particularly prominent in the middle A ring; almost 
all the the radial scans taken at the equinox show a temperatnre peak aronnd 
129,000 km for the A ring. 

In this paper, we examine two possibilities for the eqninox temperatnre 
anomaly of the A ring. The first one is that the mnlti-particle-layer structnre 
assnmed in Morishima et al. (2009) is inappropriate for the A ring. Since 
the radial peak location of the eqninox temperatnre coincides with the peak 
location of the amplitnde of photometric azimnthal brightness asymmetry 
caused by self-gravity wakes (Dones et ah, 1993; French et ah, 2007), the 
discrepancy between the modeled and observed eqninox temperatnres may 
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be due to wakes that are not taken into account in Morishima et al. (2009). 
In Morishima et al. (2014), we have developed a new thermal model in which 
wakes are represented by inhnitely long elliptical cylinders, originally intro¬ 
duced by Hedman et al. (2007). This model can reproduce the observed 
azimuthal temperature modulation, which is found to be caused by the vari¬ 
ation of the geometric hlling factor seen from the Sun. In the present paper, 
we apply the wake model of Morishima et al. (2014) to the A ring data 
at the equinox. We also perform some additional parameter studies using 
the multi-particle-layer model of Morishima et al. (2009), because the entire 
parameter space has not been exhausted in Spilker et al. (2013). 

The second possibility, for which we spend much more time than the hrst 
one, is incomplete cooling down of A ring particles at the equinox due to 
the effect of a seasonal thermal inertia. If the thermal inertia of the interior 
of a ring particle is high, the interior cools more slowly than the particle 
surface. The heat transport from the interior warms up the surface, which 
leads to a higher ring temperatnre than a model prediction withont this 
seasonal effect. Therefore, if the eqninox temperature anomaly is really due 
to the seasonal effect, we are able to give constraints on the seasonal thermal 
inertia, using seasonal temperature variation curves including the equinox 
data. The thermal skin depth associated with the seasonal effect is a few 
order of magnitude larger than the dinrnal skin depth, and is comparable to 
the ring particle size. Thus, combining with the diurnal thermal inertia, it is 
possible to clarify how the flufhness of ring particles changes with their depth. 
The seasonal thermal inertia can be constrained only if the particle size is 
larger than the seasonal thermal skin depth. Otherwise, the temperatnre 
of the entire particle down to its center qnickly adjusts to the equilibrium 
temperature determined from external heat sources. Thus, either the particle 
size or the seasonal thermal inertia would be constrained. 

In Section 2, we explain our methodology. In Section 3, we analyze the 
CIRS eqninox data, using two sophisticated radiative transfer models which 
assnme extremely different ring structures. We will show that the equinox 
temperatnres for the A ring are much higher than the model estimates re¬ 
gardless of ring structure assumed in the models, as long as only the flnx from 
Satnrn is taken into account. In Section 4, we analyze seasonal temperature 
variation using all the existing CIRS data. We introduce a simple seasonal 
temperatnre model in which the thermal diffusion equation is nnmerically 
solved for the temperature evolution of the particle interior. We adopt two 
models for the internal structure of a particle: the hrst one assnmes that 
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Fig. 1. Schematic diagram of the methodology. 


the internal density and the thermal inertia are uniform with depth whereas 
the second one assumes that a particle has a dense core covered by a fluffy 
regolith mantle. The particle size, the seasonal thermal inertia, and the 
core size estimated from model hts are presented. In Section 5, we discuss 
implication from our results. The conclusions are given in Section 6. 

2. Methodology 

We briefly review our methodology (Fig. 1) to familiarize the reader with 
the highlights from previous papers that are applicable to this work. In the 
present paper, we analyze two different data sets, the equinox data and the 
seasonal data, using three different thermal models: the multi-particle-layer 
model (Morishima et al., 2009), the mono-layer wake model (Morishima et 
al., 2014), and the seasonal model. 

As mentioned in Section 1, the hrst two models that precisely conserve 
the energy are used for comparison with the equinox data (Section 3). Since 


6 






both models explicitly give the temperature distribution of the particle sur¬ 
face (or the wake surface) by solving the radiative heat transfer problem, 
they are applicable to any observational geometries. These models also take 
into account diurnal temperature variations by solving the thermal diffusion 
equation for the thin layers near the particle (wake) surface. However, these 
models calculate energy balance of the ring at a hxed solar elevation angle, 
neglecting effects due to the seasonal thermal inertia. Including the seasonal 
effects in these two models is beyond the scope of the present paper, partly 
because of difficulties in coding. Even if we develop such models, simulations 
are most likely computationally very intense because we need to perform 
radiative transfer calculations at every solar elevation angle. 

Instead, we introduce a simple seasonal model for analysis of the seasonal 
data (Section 4). The model ignores diurnal variation of the incoming flux 
but takes into account the seasonal effect by solving the thermal diffusion 
equation for the entire particle down to its center. The model calculates the 
spherically averaged temperatures of a ring particle without considering a 
temperature distribution over the particle surface. The spherically averaged 
solar flux onto the particle is given by the analytic formula in Altobelli et 
al. (2008). The formula takes into account mutual shading assuming a 
mono-particle-layer ring. In the model, mutual heating is also represented 
by a simple analytic formula (Froidevaux, 1981; Ferrari et ah, 2005). Use 
of these analytic formulas allows us to omit any computationally intense 
radiative transfer calculations. The problem is that the simple seasonal model 
itself does not give an accurate equinox temperature with the Saturn flux 
only, while it is essential in the estimation of the seasonal thermal inertia. 
Therefore, we use the equinox temperatures derived by the two radiative 
transfer models (Morishima et ah, 2009, 2014). Fortunately, as we will show, 
both radiative transfer models give similar equinox temperatures even though 
very different ring structures are assumed. Thus, the equinox temperature 
for the mono-particle-layer model is expected to be similar as well. 

Readers familiar with these previous studies may want to skip Section 3 
and start reading from Section 4. 
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3. Comparison between modeled and observed equinox tempera¬ 
tures 

3.1. Equinox data 

The solar equinox occurred on August 11th in 2009. During the equinox 
period (August llth-13th), CIRS obtained 15 radial scans (Spilker et ah, 
2013). The absolute solar elevation angle dnring the eqninox period was 
between 0 and 0.036°. The very low solar elevation angle guarantees that 
the direct solar illnmination is negligible. The spacecraft elevation angles 
were about 20° for all the scans. 

We analyze the equinox data in two forms. The first one is mean temper- 
atnres observed aronnd 129,000 km (±500 km) with various observational 
geometries. The equinox temperature at a certain distance from Saturn in¬ 
creases with decreasing Saturn phase angle (Saturn-ring-observer angle) in¬ 
dicating that ring particles are not isotropic emitters, althongh the geometry 
dependence is very weak for the A ring. Instead of the Satnrn phase angle, 
we use the spacecraft hour angle in the rotating frame, Ts/c,R; introdnced in 
Spilker et al. (2013), to examine the temperature asymmetry between the 
leading and trailing sides of the ring. The angle is give as Ts/c,r = -^s/c~-^R; 
where Ts/c is spacecraft hour angle around a footprint and Lr is the ring 
local honr angle around Saturn. In this frame, the direction of Saturn seen 
from the ring is always at Ts/c,r = 12 h (= 180°), and the leading and 
trailing sides of the ring are seen when Ts/c,r < 12 h and Ts/c,r > 12 h, 
respectively. Since the spacecraft elevation angles are almost the same for all 
the equinox scans, the Saturn phase angle simply increases as Ts/c,r recedes 
away from 12 h. The second one is the temperatures at various saturnocen- 
tric distances. We choose two radial scans that give the highest and lowest 
mean temperatures at the equinox as representatives. Radial temperatnre 
prohles of other scans are very similar to these two scans. More details for 
the eqninox data set are found in Spilker et al. (2013). 

3.2. Radiative transfer models 

We test whether thermal models can reprodnce the equinox tempera¬ 
tures. Two radiative transfer models are employed: a multi-particle-layer 
model (Morishima et ah, 2009) and a mono-layer wake model (Morishima et 
ah, 2014). In the models, Saturn’s thermal emission and solar illumination 
reflected by Satnrn are taken into account, and direct solar illumination is 



ignored. The effective surface temperature of Saturn is 96.67 K (Li et ah, 
2010) and Saturn’s oblateness is 0.1. 

3.2.1. Multi-particle-layer model 

In the multi-particle-layer model of Morishima (2009), the vertical thick¬ 
ness of a ring is assumed to be much larger than the particle size. The model 
numerically solves the classical radiative transfer equation (Chandrasekhar, 
1960; Modest, 2003) using the multi-stream method (we use 60 streams). 
We adopt the plane-parallel approximation so rings are horizontally homoge¬ 
neous. The model takes into account heat transport due to particle motion in 
the vertical and azimuth directions by following temporal thermal evolution 
of thousands of ring particles. It assumes, instead of an actual continuous 
size distribution, a bimodal size distribution consisting of fast and slow ro¬ 
tators. Fast rotators are small-rapidly spinning particles having spherically 
symmetric thermal structure whereas slow rotators are large non-spinning 
Lambertian particles. The ring effective temperature is derived by summing 
up the thermal emission of all ring particles toward the observer, taking into 
account the attenuation factor due to other particles on the line of sight. 

In the present work, we consider two extreme cases of particle spins: fast 
rotators only and slow rotators only. In the model, we adopt the following 
parameters: the infrared emissivity cir = 1.0, the thermal inertia T = 10 
Jm“^K“^s“^/^, and the bolometric Bond albedo = 0.5. The latter two 
parameters are estimated by Morishima et al. (2010, 2011). 

3.2.2. Mono-layer wake model 

In the mono-layer wake model of Morishima et al. (2014), the temperature 
distribution over the wake surface is calculated instead of considering the 
temperatures of individual particles. Wakes are approximately represented 
by inhnite, parallel, regularly spaced, opaque elliptical cylinders, with vertical 
thickness FT, horizontal width hF, and wavelength A, as adopted in Hedman 
et al. (2007). Our model ignores interwake particles for simplicity. The 
surface of an elliptical cylinder is divided into multiple facets and interactions 
between facets (mutual heating and multiple scattering of photometric light) 
are calculated using a radiosity method. As discussed in Morishima et al. 
(2014), it is likely that some heat transport occurs within wakes; this may 
be due to vertical motion of particles, or particle spins may play a role if 
wakes are represented by a mono-particle-layer. To represent heat transport 
approximately, the incoming flux is smoothed over the wake surface and 
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A ring; a = 129,000 km 



0 5 10 15 20 

LHAg (Hours) 


Fig. 2. Equinox temperature variation with local hour angle of Cassini in the 
rotating frame. The saturnocentric distance of the ring is 129, 000 km and 
the radially averaged (±500 km) temperatures are shown. The direction of 
Saturn is 12 hr in this frame. Black pluses (with error bars) are the observed 
temperatures and blue (fast rotators) and red (slow rotators) squares are 
temperatures from the multi-particle-layer model. 

its degree /mix (0 < /mix < 1) is given as a parameter. If /mix = 1, the 
wake surface is isothermal. The ring effective temperature is calculated by 
summing up the thermal emission of facets seen from the observer. 

We use /mix = 0.8, the wake bolometric reflectance Ay = 0.35, and T = 10 
Jm“^K“^s“^/^ as estimated by Morishima et al. (2014) unless otherwise 
noted (/mix is varied in Fig. 4). The pitch angle of wakes is assumed to be 
30°. A similar value is derived by Ferrari et al. (2009). 

3.3. Equinox temperature results 
3.3.1. Munlti-particle-layer model results 

Figure 2 shows the variation of the equinox temperatures with Ts/c,r at 
a = 129, 000 km. The observed temperatures in the hgure are the same with 
those shown in Spilker et al. (2013) (their Figs. 3 and 7). In the model, 
we adopt the ring optical depth r of 0.5. Both cases of particle spins give 
equinox temperatures that are much lower than the observed temperatures 
although slow rotators give slightly higher temperatures than fast rotators. 
Similar results were previously presented in Spilker et al. (2013) in which 
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Fig. 3. Radial profiles of equinox temperatures. Individual diamonds are 
temperatures retrieved by individual spectra. Two scans that give the highest 
(red diamonds) and lowest (bine diamonds) mean temperatnres are shown. 
Radial temperatnre prohles calculated by the multi-particle-layer model are 
also plotted by solid (slow rotators) and dashed (fast rotators) lines. 

the total cross sections of fast and slow rotators are assnmed to be ronghly 
the same, as snggested by Morishima et ah (2010, 2011). 

Fignre 3 shows the radial variation of the eqninox temperatnres at the two 
different geometries. In the model, we nse the optical depths from Voyager 
PPS (Esposito et al., 1983) smoothed over radial bins with a 1,000 km width, 
avoiding the Encke division. The observed radial prohles show a temperatnre 
peak aronnd 129,000 km whereas the model predicts that the temperatnre 
gradnally decreases with increasing distance. The ring normal optical depth 
is known to depend on observational geometry due to wakes (Colwell et ah, 
2006, 2010; Hedman et ah, 2007). If we use larger optical depths given by 
some of Cassini occultations, slightly lower temperatnres are obtained dne 
to an enhanced mntnal shading effect while the relative radial temperatnre 
variation is almost nnchanged. 

Overall, the model cannot reprodnce not only the absolnte temperatnre 
bnt also the radial prohle of the eqninox temperatures. If we adopt an in¬ 
frared emissivity lower than nnity, that can probably give slightly higher 
temperatures. Since the current mnlti-particle-layer code does not include 
scattered thermal light, we cannot directly check cases with cir < 1.0. On 
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A ring: a = 129,000 km 
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Fig. 4. Same as Fig. 2, but temperatures from the mono-layer wake model 
are compared with the observed temperatures. Dependences on the thermal 
emissivity cir, and the mixing efficiency /mix are examined. The height-to- 
width ratio, H/W, is 0.4. 


the other hand, we can roughly infer its effect from the wake model which 
includes scattered thermal light. We see that the equinox temperatures for 
riR = 0.9 are higher than those for cir = 1.0 by 1-3 K (Fig. 4). This difference 
is still insufficient to explain the observed temperature anomaly, particularly 
in the middle A ring. 

3.3.2. Mono-layer wake model results 

Figure 4 shows results of the wake model on top of the observed data 
similar to what is shown in Fig. 2. In the model, we adopt WjX = 0.4 
and H/W = 0.4. The wake width gives the normal ring optical depth of 
T = — log(l — W/X) = 0.51 which is close to that for the A ring and H/W = 
0.4 is the value suggested by Ferrari et ah (2009). In Fig. 4, we examine 
dependence on eiR and /mix- If cir = /mix = 1, the wake temperature becomes 
independent of observational geometry at a certain local hour angle. For 
/mix < 1, the observed ring temperature increases towards Ts/c,r = 12 h, 
because the surface portions facing towards Saturn have relatively higher 
temperatures. However, the temperature averaged over Ts/c,r is insensitive 
to /mix- With decreasing ejR, the ring physical temperature increases due 
to inefficient thermal emission and a part of infrared emission from Saturn 
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A ring: a = 129,000 km 
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Fig. 5. Same as Fig. 4, but the dependence on height-to-width ratio H/W is 
examined in the model. We adopt eiR = 0.9 and /mix = 0.8. 



Fig. 6. Same as Fig. 3, but temperatures from the mono-layer wake model 
are compared with the observed temperatures. The results for H/W = 0.2 
(solid) and 0.5 (dahsed) are shown. We adopt cir = 0.9 and /mix = 0.8. 
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is scattered by wakes toward the observer. In addition, a low cir causes 
a strong observational-geometry dependence of the temperature because of 
scattered thermal light. We hnd that to reconcile with the observed geometry 
dependence, cir needs to higher than ~ 0.9. The high emissivity is also 
indicated by spectroscopy (Morishima et ah, 2012), and is due to strong 
absorption of water ice at far-infrared wavelengths. As long as cjr > 0.9, 
the wake model gives equinox temperatures that are much lower than the 
observed temperatures. 

Although H/W ~ 0.4 is favored for the elliptical cylinder wake model 
and CIRS data in the middle A ring (Ferrari et ah, 2009; Morishima et ah, 
2014), estimated values of H/W vary with different observations and wake 
shape models (e.g., Cowell et ah, 2006; Dunn et ah, 2007; French et ah, 2007; 
Hedman et ah, 2007). Thus, we check how the equinox temperature varies 
with H/W (Fig. 5). It is found that the equinox temperature increases with 
H/W as long as H/W < 0.5. For a larger H/W, little increase is seen. The 
effective wake cross section seen from Saturn increases with increasing H/W, 
as long as other nearby wakes do not hinder illumination from Saturn. That 
is the case for H/W < 0.5. For H/W > 0.5, mutual shading starts to be 
effective. Similar results are found for cases in which solar illumination is the 
dominant heat source (Morishima et ah, 2014). 

Figure 6 shows radial prohles of the equinox temperature derived from 
the wake models for H/W = 0.2 and 0.5 plotted on top of the observed data 
similar to what is shown in Fig. 3. Again, the optical depth prohle from 
Voyager PPS is used. It is found that the mono-layer wake model also gives 
much flatter radial prohles of the equinox temperatures than the observed 
ones. Colwell et al. (2006) showed that H/W increases from ~ 0.2 at the 
inner A ring to ~ 0.5 at the outer A ring. If we adopt such a radial variation 
of H/W, the equinox temperature becomes almost constant over the A ring, 
inconsistent with the observations. 

We conclude that neither the wake model nor the multi-particle-layer 
model can reproduce the observed high equinox temperatures and the radial 
temperature prohles of the A ring with a plausible range of parameters. 

4. Particle properties inferred from seasonal temperature varia¬ 
tions 

In Section 3, we showed evidence of incomplete cooling down of Saturn’s 
A ring at the equinox. This fact allows us to give constraints on particle 
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properties such as the size and the seasonal thermal inertia. In this section, 
we develop a simple time-dependent model for seasonal ring temperature 
variations and compare model results with observations. 

4 . 1 . Data selection 

As of the middle 2014, we have about 2,000 CIRS scans. Each scan usnally 
contains thousands of spectra (and temperatnres). We make data selections 
from the entire set of scans as follows. We choose data with moderately good 
spatial resolntions to examine radial variation of particle properties. Specif¬ 
ically, data with a radial length of a footprint projected on the ring plane 
of less than 5,000 km are selected. While our seasonal model does not take 
into acconnt observational geometry dependence, the data show variations of 
temperatnre with geometry, particnlarly when \B'\ is high. If \B'\ < 2°, we 
do not remove any data as the geometry and dinrnal dependences of temper¬ 
atnres are small enongh. If \B'\ > 2°, we remove data at solar phase angles 
lower than 50° as the low phase data show high temperatures. It is more 
appropriate to calculate the mean temperatnre by integrating the thermal 
flux over the solid angle as done in Pilorz et al. (2014) for the B ring. Never¬ 
theless, since the A ring shows relatively weak phase dependences (Altobelli 
et ah, 2009), this simple removal of low phase data has only a minor effect. 
For \B'\ > 2°, we also remove data near Satnrn’s shadow as the temperatnres 
in the shadow are mnch lower than those aronnd the noon. We use data only 
with 100° < Lr < 300°. 

The temperatures in a single scan (nsnally a radial scan) are averaged over 
14 radial bins with a width of 1,000 km and the bin centers from 123,000 
km to 136,000 km. The nnmber of data points in a radial bin for a single 
scan is generally about 10. The averaging is done after removing the data 
ontside 3a, where a is the standard deviation of the temperatnre. If a after 
removing noisy data is still larger than 3 K, we do not use the scan for this 
radial bin. 

We make further data cleaning. Using the model introdnced in Sec¬ 
tion 4.3, we derive the best-fit curves for seasonal temperatnre variations. 
If the difference between the best £t temperatnre and the mean observed 
temperatnre is larger than 5 K, we remove this scan. These abnormal data 
most likely snffer from some sort of instrnmental errors. Since these data are 
rather rare, model fits are not altered regardless of them. 

One may find that we have a great deal of CIRS data before the eqninox 
whereas mnch less data are available after the eqninox (e.g.. Fig. 7). After 
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Temperature vs. Solar elevation angle 



Solor elevation angle (deg) 


Fig. 7. Seasonal temperature variation around the solar equinox {B' = 0) 
at 129, 000 km. The blue squares and red diamonds are the southern and 
northern face temperatures, respectively. The solid line is the time since the 
equinox. The temperatures in between dashed and dotted vertical lines are 
used for measurements of the asymmetry in Fig. 8. 


the equinox, Cassini changed its orbital inclination relative to the Saturn’s 
equatorial plane from moderately high values to nearly zero to observe Sat¬ 
urnian satellites. Only a small number of ring observations were made soon 
after the equinox, as observational geometries were not suitable for the main 
rings. The lack of post-equinox data unfortunately causes large uncertainties 
in parameters estimated in the present work. 

4-2. Asymmetry in seasonal temperature curves around equinox 

Before moving to modeling, we will show additional evidence of incom¬ 
plete cooling down at the equinox, using seasonal data only. Incomplete 
cooling down means that it takes hnite time for ring particles to cool down 
due to a hnite thermal inertia and a hnite size. This should also mean that 
it takes hnite time for particles to be warmed up by direct solar illumination 
after the equinox. As a result, it is expected that the ring temperature before 
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Fig. 8. Temperature asymmetry around the equinox. The upper and lower 
panels show the temperatures before the equinox subtracted by those after 
the equinox in the ranges of \B'\ = 2.1-2.5° and 0.15-0.35°, respectively. The 
blue and red marks are for the lit and unlit faces, respectively. 


the equinox is higher than that after the equinox for the same absolute solar 
elevation angle. This asymmetry in seasonal temperature curves is indeed 
identified. 

Figure 7 shows the seasonal temperature variations of the middle A ring 
(at 129,000 km) over ~ a year around the equinox. The blue symbols are 
temperatures of the southern face of the ring the red symbols are tempera¬ 
tures of the northern face. We call a face directly illuminated by the Sun the 
lit face and the other face the unlit face. The southern face changes from the 
lit face to the unlit face at the equinox. Asymmetries for both the lit and 
unlit faces are examined (e.g., comparison between the pre-equinox southern 
face and the post-equinox northern face). Figure 7 clearly shows that the 
pre-equinox temperatures are larger than the post-equinox temperatures at 
the same absolute solar elevation angle both for the lit and unlit faces. 

To quantify asymmetry of temperatures, we take averages of temperatures 
in the intervals between \B'\ = 2.1° and 2.5° and between \B'\ = 0.15° 
and 0.35°, and subtract the mean temperature after the equinox from the 
mean temperature before the equinox for the same absolute solar elevation 
angle. Figure 8 shows the temperature differences at various saturnocentric 
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distances for the intervals between \B'\ = 2.1° and 2.5° (the upper panel) 
and between \B'\ = 0.15° and 0.35° (the lower panel). The blue symbols are 
for the unlit face and the red symbols are for the lit face. It is found that 
the asymmetry is particularly large at the middle A ring, where the observed 
equinox temperatures are much larger than the modeled temperatures. The 
large temperature asymmetry, in addition to the large anomaly in the equinox 
temperature (Section 3), strongly indicates the efficiency of cooling down is 
low in the middle A ring. On the other hand, the asymmetry is relatively 
small at the inner and outer parts of the A ring, where relatively small 
anomalies in the equinox temperatures are seen. Thus, it is likely that some 
particles properties vary radially across the A ring. Those will be addressed 
in the following sections. 

4-3. Seasonal model 

4.3.1. Energy balance for ring particles 

We introduce a seasonal temperature variation model. The model takes 
into account time-dependence by solving the thermal diffusion equation of a 
particle. In the model, the ring temperature is represented by the surface 
temperature of a particle and the ring particle temperature is assumed to 
be spherically symmetric. Therefore, dependency of ring temperature on ob¬ 
servational geometry is not considered although the southern and northern 
face temperatures are calculated separately by giving different solar fluxes to 
these faces. The diurnal variation of the incoming flux is also ignored. We 
adopt the analytic formula of the spherically averaged direct solar flux given 
by Altobelli et al. (2008). This formula takes into account mutual shading 
of particles and accurately reproduces the flux onto a mono-particle-layer 
ring (Froidevaux, 1981). The model of Froidevaux (1981) is known to accu¬ 
rately reproduce the temperature variation with solar elevation for all rings 
at any observational geometries with an appropriate choice of the geometry- 
dependent particle albedo (Flandes et ah, 2010). The model requires an 
accurate Saturn flux that determines the equilibrium equinox temperature 
that is obtained when the time-dependence is ignored. For an accurate Sat¬ 
urn flux, we use the equinox temperatures given by our own radiative transfer 
models, instead of the approximate formula given by Froidevaux (1981). We 
use the equinox temperatures averaged over the four radial prohles shown in 
Fig. 3. These temperatures are calculated by the multi-particle-layer model 
but the mono-layer wake model gives similar temperatures (Fig. 6). 
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To develop a model of the seasonal temperature variation of a ring, we 
first need to determine the heliocentric distance and the solar elevation 
angle B' as a function of time, t. We simply assume that Saturn has a 
Keplerian orbit around the Sun. The details are given in Appendix A. 

Time evolution of the particle temperature T is solved by the one-dimensional 
thermal diffusion equation: 

A . J_iT , (1) 

dt pc C dr \ dr ) ' 

where K is the thermal conductivity, p is the internal density, and C is the 
specihc heat. We use a hxed value, C = 450 J kg“^ K“^, which is the value 
of water ice at ~ 50 K (Shulman, 2004), whereas K and p are varied in the 
model. We ignore the temperature dependence of C and K for simplicity. 
This is probably a safe assumption as the time-dependence is important only 
around the equinox and the temperature variation during that period is not 
large. At high solar elevation, the ring temperature is very close to the 
equilibrium temperature given with A' = 0. 

The boundary condition at the surface is given by 

^ -^Sun + Tsat - eiRcrT > (2) 

where Fsun and Fsat are the fluxes from the Sun and Saturn, respectively, ctsb 
is the Stefan-Boltzmann constant, and flp is the solid angle of surrounding 
particles. The third term in the right hand side of Eq. ([2]) is the radiation 
flux from the particle, and a part of the emission is blocked by surrounding 
particles that are assumed to have the same temperature as the subject 
particle. We adopt ffp = 6[1 — exp(—r)] (Ferrari et ah, 2005), where r is the 
ring optical depth. 

The solar flux is given by 

(l^) ■f’®' (3) 

where A is the albedo of the particle (either or A „niH , dehned below), / 
is the spin factor, Fe = 1370 Wm ^ is the solar flux at lAU. We introduce 
two albedos: the lit-face albedo A^t that is used when the face is directly 
illuminated by the Sun and the another one for the unlit face Auniit- It is 
expected that Am is close to the bolometric Bond albedo of the particle. On 
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the other hand, Amii t, includes the effect due to vertical heat transport that 
decreases with increasing r. Thus, it is expected that Auniit is larger than 
A\it as also shown by Flandes et ah (2010). The model assumes that the 
timescale of vertical heat transport is short compared with the orbital period 
of Saturn, Ts. This is a reasonable assumption if vertical heat transport is 
taking place by particle motion (Morishima et al., 2010; Pilorz et ah, 2014). 
The spin factor / is 2 for slow rotators and 4 for fast rotators. We adopt 
/ = 3 that is suggested from equinox modeling (Spilker et al., 2013). In 
Eq. ([3]), S represents the factor due to shadowing of surrounding particles 
and we use the approximate formula given by Altobelli et al. (2008): 


SiB^r) 


sin \B'\ 


1 — exp (—r) 


exp 


sin \B' 


( 4 ) 


The Saturn flux is given by 


Esat — eiRd Tj 


eq 



( 5 ) 


where Tgq is the equilibrium equinox temperature due to the Saturn flux only 
(note that T = Tgq if Fsun = K = 0). As stated above, we take Tgq from 
our radiative transfer calculations shown in Section 3. We fix cir = 1 in the 
seasonal model. Note that cir is not a free parameter as scattered thermal 
light is not taken into account in the model; dividing Eq. ([2]) by cjr, the 
meaningful parameters are found to be iP/eiR and (1 — A)/(/eiR), but not 
eiR itself. 

The time in the thermal diffusion equation can be normalized by the 
inverse of the mean motion of Saturn, The length can be normalized 

by the seasonal thermal skin depth i given by 


K 


pCujs pCy/uJ^ 


~ 0.60 


10 Jm“^K“^s“^Ay ^450 kgm 


P 


1-3 


m, 


( 6 ) 


where T = yjpCK is the thermal inertia. This length is a few order of mag¬ 
nitude larger than the diurnal skin depth, ~ 1 mm (Morishima et ah, 2009), 
and comparable to the ring particle size suggested from occultations (French 
and Nicholson, 2000; Cuzzi et ah, 2009). After the normalization of Eq. (IT]) . 
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r appears as a single parameter rather than p, C, and K, independently. If 
£ 3> i?, the thermal relaxation time (or the cooling time) of the particle is 
given by 

{47r/3)R^pCT RpC 
~ ATiRhmaT^ ~ 3eaT^ 

" (l^) (s^) <"> 

Thus, the equinox temperature sensitively depends on RpC. On the other 
hand, if i R, the thermal relaxation time is given by 

AnR^ipCT _ r 
“ AnR^ei^aT^ ~ eaT^^ 

- ^ 1 / 2 ) (®) 

Vio jm- 2 K-is-h 2 y \^5o k; 

Thus, the equinox temperature sensitively depends on T. If i? ~ £, the 
thermal relaxation time takes an intermediate value between Treii and Trei 2 
and the equinox temperature depends on both RpC and T. 

4-3.2. Particle internal structure models 

For the interior of the particle, we consider two models (Fig. 9). The hrst 
model, called the uniform model, adopts uniform F and p with depth. We 
adopt p = 450 kgm“^. This value is suggested from photometric azimuthal 
brightness asymmetry (Salo et al., 2004, French et ah, 2007). The parameters 
to be htted are F, R, Ant, and v4unlit- For simplicity, we consider only single¬ 
sized particles whereas actual ring particles have a size distribution (Cuzzi 
et al., 2009). The inner boundary for Eq. ([T]) is set to be a larger one of 0 
and R — 7i, and the boundary condition there is given as 


dT 



0 . 


(9) 


In the second model, called the core model, we consider a dense core 
with a radius Rc covered by a fluffy regolith mantle with a thickness of 
L = [R — Rc). We assume that the core is made of ice Ih with no porosity 
and adopt the properties at 50 K: the internal density pc = 933.74 kgm“^ 
(Feistel and Wagner, 2006) and the thermal conductivity = 12 W m“^ 
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p = 450 kg m ^ 
r = const 


Uniform model 



Core model 


Core 

p = 933.74 kg m ^ 
r = 2245Jm2Kisi/2 


Regolith mantle 
p = 100 kg m'^ 
r = 10 Jm 


Fig. 9. Schematic illustration of two internal structure models. 
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(Andersson and Inaba, 2005). These values give the thermal inertia of the 
core as T^ = 2245 The regolith mantle is assumed to have p = 

100 kgm“^ and T = 10 Jm“^K“^s“^/^; the value of T is those estimated from 
eclipse cooling (Morishima et ah 2011, 2014) and p is those indicated from the 
relationship between the porosity and the thermal conductivity (see Fig. 14 
for details). In the core model, R and L/R are varied as htting parameters. 
We explore only in the range of 0.05 < L/R < 0.5, as a large(small) L/R 
gives a mean internal density that is too low(high). We use A\\t and A „nii t 
from the uniform model. Since the core’s K is large enough, the thermal 
diffusion time scale for a meter-size core is short enough (less than a day). 
Thus, the spatial temperature variation inside the core can be ignored and 
the core can be well represented by a single grid in numerical calculations. 
The equation of the core temperature change is given by 


R^pC dT 
3 dt 



1 

r=Rc 


( 10 ) 


where the thermal conductivity, at the core-mantle boundary is set to 
be 2{K~^ + ~ 2K (Shoshany et ah, 2002). 

For the uniform model, the particle or the surface portion is divided 
into 100 layers with an identical thickness, and the temperature change is 
numerically calculated using the Crank-Nicholson method (Press et a/., 1986; 
Chap. 19.2). After convergence tests, we adopt the timestep of 10“®Ts. This 
very small timestep is required because the solar flux changes very rapidly on 
a short timescale around the equinox. For the core model, the regolith mantle 
is divided into multiple layers. It is found that numerical instabilities occur 
if the thickness of each layer is too small as compared with the skin depth. 
After some experiments, we adopt the number of layers, N-^, as follows: 

TVl = min(2200y,100). (11) 


In the seasonal model, we again use the optical depths from Voyager PPS 
(Esposito et ah, 1983) smoothed over each radial bin, avoiding the Encke 
division. We also used some of prohles taken by Cassini UVIS (Colwell et ah, 
2010), but we did not see any signihcant differences in the htted parameters. 
Thus, those results will not be presented. 


4.4- Fitting procedures 

For model htting, we divide the observed data into four portions: (1) the 
southern face data before the equinox, (2) the southern face data after the 
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equinox, (3) the northern face data before the equinox, and (4) the northern 
face data after the equinox. The total number of scans for the portion (1) 
is dehned to be nsi,max, and the numbers for other portions are dehned as 
well. The Usi-th scan in the portion (1) contains nfi,max(nsi) temperatures 
in a designated radial bin of the width of 1,000 km. The mean temperature 
and the standard deviation of the Ugi-th scan are dehned to be Tobsi(nsi) and 
(Jobsi(nsi). The reduced chi-square for the portion (1), Xi, is then calculated 
as 


Xi(r , R, y4iit, Aunlit) — 


nsl,r 


-M 


X 




71=1=1 


^obsl (nsi) ^siml (nsi) T, Rj idmjiit) 
l^obsl (nsi) 


( 12 ) 


where M = 4 is the number of the parameters to be htted, Tgimi is the 
modeled temperature, and w(nsi) is the weight given to each scan. After 
some experiments, we determined to adopt the following form for the weight: 


w(nsi) = CiUfi 

,max (usi) exp 


\B'\ 
3 deg 


(13) 


where B' is the solar elevation angle at the time when the scan was taken 
and the coefficient ci is derived by normalization of the weight as 


^sl ,max 

w{nsi) = 1. 

71sl = l 


(14) 


Eq. flT^ means that more weight is given to data near the equinox. The 
temperature at high \B'\ is almost the equilibrium temperature obtained by 
ignoring the time dependence, and is primarily adjusted by the albedos A\\t 
and A ,,nib ,. Therefore, we cannot well constrain T and R if we weight all the 
data equally. 

The reduced chi-squares for the portions (2)-(4) are calculated in a sim¬ 
ilar manner. Then, the reduced chi-square for the uniform model is simply 
dehned by averaging the four values as 

xi = ^ (xi + X 2 + Xs + xl) ■ (15) 

A set of the best-ht parameters are obtained at the chi-square minimum, 
x\ min, in the four-dimensional parameter space. The error bars are estimated 
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from the maximum and minimum values of the parameters on the contour 
of xi/XAmin = 2. A similar procedure gives the reduced chi-square, Xm for 
the core model and the best-£t parameters are obtained at the chi-square 
minimum, 

4-5. Uniform model results 

Before discussing the best-fit parameters, we show how seasonal temper¬ 
ature curves change with different input parameters. Figure 10 shows the 
observed and modeled temperatures near the equinox at 129,000 km. In the 
model, the particle size and the thermal inertia are varied while the best-£t 
values of the albedos. Am = 0.498 and Am = 0.625, from Fig. 11 are used. 
The upper panel shows the case of F = 10 Jm“^K“^s“^/^. The seasonal ther¬ 
mal skin depth i is 0.60 m in this case. The seasonal temperature variation 
is almost independent of R as long as R > i. For the case of R = 0.3 m, the 
equinox temperature is slightly lower than large R cases. Both the lit and 
unlit face temperatures are nearly symmetric with respect to the equinox, 
as the thermal relaxation time, rrei 2 [Eq. ([H])] (= 198 days), is short enough. 
The timescale ri.ei 2 roughly coincides with the duration of \B'\ < 1°, and 
it is expected that the ring temperature is nearly the equilibrium tempera¬ 
ture for \B'\ > 1° in this case. The lower panel shows the case of F = 50 
Jm“^K“^s“^'^^ and this case gives i = 3.0 m. In this case, the equinox tem¬ 
perature depends on R, as R < 1. The lit and unlit face temperatures are 
asymmetric with respect to the equinox particularly for R = 1.0 and 3.0 m; 
the temperature after the equinox is lower than that before the equinox at 
the same \B'\. The thermal relaxation time, rrei 2 is 993 days, and this is 
roughly the duration of |i?'| <6 — 7°. Indeed, we see convergence of the tem¬ 
perature curves for three different values of R when \B'\ > 6 — 7° (outside 
of the hgure). For R = 0.3 m, the temperature curves show only a small 
dependence on F and this indicates that the thermal relaxation time is given 
by Treii (33 days). The modeled equinox temperatures are lower than the 
observed ones if F = 10 Jm“^K“^s“^/^ or R = 0.3 m. Thus, F and R are 
expected to be larger than those values at 129,000 km. 

Figure 11 shows the best-£t values of F, R, Am, and Auniit- The reduced 
chi-square minimum, Xmin A is also shown. The thermal inertia F is typically 
30-50 Jm“^K“^s“^A the middle A ring, while it is ~ 10 Jm“^K“^s“^A 
in the inner and outermost A ring. The latter value is close to the diurnal 
thermal inertia derived from eclipse temperatures (Morishima et ah, 2011, 
2014). The particle size is 1 -2 m in most of saturnocentric distances but 
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Solar elevotlon angle (deg) 


Temperoture vs. Solar elevation angle 



Solar elevation angle (deg) 


Fig. 10. Examples of seasonal temperature variations for the uniform model 
at 129, 000 km. The observational data are the same as those in Fig. 7. The 
upper and lower panels show the cases of F = 10 and 50 Jm“^K“^s“^/^. The 
results of three different radii R = 0.3(solid), l.O(dotted), and 3.0 m (dashed) 
are shown for the southern (blue) and northern (red) faces. 
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Fig. 11. Estimated parameters for the uniform model: (a) thermal inertia F, 
(b) particle radius R, (c) particle albedos and Amu t, and (d) goodness 
of hts Xa min- The ranges of F and R are limited a priori in data hts as 4 
Jm“^K“^s“^/^ < F < 181 Jm“^K“^s“^/^ and 0.125 m < i? < 16 m. 
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Temperature vs. Solor elevation angle 



Solor elevotion ongle (deg) 


Temperature vs. Solar elevation angle 




Fig. 12. Data fitting with the uniform model at 125, 000 km. (a) Best-ht 
temperature curves as a function of B'. The northern and southern face 
temperatures are represented by red and blue colors, respectively, (b) En¬ 
largement of the region surrounded by dotted lines in (a), (c) Contour of 
xi/ximin — 1 on the F vs. R plane sliced at the best-ht Awi and Aunut. The 
point of the best-ht F and R is shown by the plus mark. 















































Temperoture vs. Solor elevation ongle 



Temperature vs. Solar elevation angle 



Solar elevation angle (deg) 



Fig.13. Same as Fig. 12, but at a = 129,000 km. 
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larger values are seen in the inner A ring. It is found that the error bars of F 
and R are very large and the upper limits of these parameters are not usually 
well constrained except for some cases such as F at 136,000 km. The value 
of Aiit is about 0.4-0.5 (note again that we adopt the spin factor / = 3), 
consistent with those estimated by Morishima et al. (2010). The value of 
Auniit is about 0.6. The ratio of Aunut to A^t is consistent with those of 
Flandes et al. (2010). In all radial locations, ximin < 1) indicating that the 
hts are good. 

Examples of the best-£t temperature curves are shown in Figs. 12 (a = 
125,000 km) and 13 (a = 129,000 km). The fits are found to be good for the 
entire range of B’. The temperature asymmetry near the equinox is caused by 
the effect of incomplete cooling down as we discussed, while the asymmetry 
at high \B'\ is caused by different saturnocentric distances. Figures 12c and 
13c shows the contours of xi/XAmin — 1 on the plane of F vs. R sliced at 
the best-£t values of Ant and A i,nii t,. As shown in Fig. 11, F and R are not 
well constrained independently. However, the region of small Xa values is 
conhned in a narrow band on the plane of F vs. R. We can clearly exclude 
particle with both high F and R. The lower limits of F and R are also well 
constrained. The contour for a = 129,000 km has a band of low Xa which 
is shifted toward a upper-right position compared with that for a = 125,000 
km. This means that at least one of F or i? needs to be larger at a = 129,000 
km than a = 125, 000 km. French and Nicholson (2000) showed that the size 
distribution of ring particles is similar across the A ring inside the Encke gap 
(a= 133,589 km). Thus, the radial variation of R is probably small {R is 
likely to be 1-2 m even for the inner A ring), whereas F should be relatively 
large in the middle A ring compared with the inner A ring (see Section 5 for 
more discussions). 

4.5.1. Expected regolith density 

The seasonal thermal inertia estimated in our model (Fig. 11) is found 
to be comparable to or somewhat larger than the diurnal thermal inertia. 
The seasonal thermal inertia is still much lower than that for rigid water ice, 
~ 2200 Jm“^K“^s“^/^, indicating that regolith of ring particles is likely to 
be very fluffy down to their deep interiors. Many studies have been made 
of the relationship between the thermal conductivity, K, and the porosity, 
p. The exact relationship varies from author to author, but in principle K 
signihcantly decreases with increasing p. Figure 14 shows two examples that 
are likely to be relevant to Saturn’s rings; one derived from numerical sim- 
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Fig. 14. The thermal conductivity K as a function of porosity, p. The conduc¬ 
tivity is normalized by its value at p = 0. The numerical model of Shoshany 
et al. (2002) gives a fitting curve, K/K{p = 0) = (1 — 
where po = 1 ~ (1 ~ pY^^. In their model, the ratio of the largest grain size 
to the smallest grain size is given by and we adopt n = 2. Krause et 

al. (2011) conducted heat conduction experiments using silicon beads and 
found the relationship K/K{p = 0) = exp (—7.91p) (their maximum p is 
0.85). The yellow rectangle indicates the estimated ranges of K and p for 
the A ring particles. For the range of p, we take p = 270-675 kg m“^ (the 
upper and lower limits are not strict) from Fig. 23 of French et al. (2007) 
and assume that ring particles are pure water ice; p = 933.74 kg m“^ at p = 
0. The range of K is from this study. 
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ulations of Shoshany et al. (2002) and another one derived from laboratory 
experiments of Krause et al. (2011). The range of K{p)/K{{)) from our model 
results is about 10“^ - 10“^ and is shown by the arrow in Fig. 14, where we 
use the values of p and C assumed in our model in the conversion from F to 
K. Although the two curves of K differ from each other, both indicate that a 
high porosity {p > 0.8) is necessary to give K as low as those estimated from 
our model. If ring particles mostly consist of water ice (Cuzzi et ah, 2009), 
this means that p ~ 200 kg m“^ or even less. Probably, ring particle regolith 
is as fluffy as fresh snow on the Earth (p ~ 100 kg m“^; Judson and Doesken, 
2000), as the regolith is likely to be impact ejecta that gently accumulated 
on particles (Cuzzi and Estrada, 1998; Elliott and Esposito, 2011, Tiscareno 
et ah, 2013; Schmidt and Tiscareno, 2013). 

However, it is known that p of A ring particles should not be that low, 
because self-gravity wakes are observed. French et al. (2007) compared the 
azimuthal brightness asymmetry observed by the Hubble Space Telescope 
with those from their iV-body ray-tracing simulations and found that the 
observed amplitude of azimuthal asymmetry is reasonably reproduced if p ~ 
450 kg m“^ (the value we assumed in the uniform model). Thus, the uniform 
model is not self-consistent as the estimated K is too low for the suggested 
value of p or p (Fig. 14). To resolve this issue, in the next section, we 
introduce the core model that allows moderately high mean internal density 
and very fluffy surface regolith. 

4 . 6 . Core model results 

We hrst examine basic trends of seasonal temperature variations with 
different input parameters, as done for the uniform model. Figure 15 shows 
the seasonal temperature curves at a = 129,000 km for L/R = 0.05 and 
0.2. Three different particle radii, R = 0.3, 1, and 3 m, are considered. The 
results may be compared with those in Fig. 10 for the uniform model. If the 
regolith thickness is large enough, the temperature curve is similar to that for 
the uniform model with F = 10 Jm“^K“^s“^/^. The case oi L/R = 0.2 (upper 
panel) is close to this limit, although the equinox temperatures are slightly 
increased due to the effect of the core. Since the core cools down slowly and is 
much warmer than the regolith mantle at the equinox, the core warms up the 
regolith mantle. With decreasing L/R, the equinox temperature increases. 
This effect is similar to increasing F in the uniform model. However, the 
shapes of the temperature curves are slightly different between the core model 
with L/R = 0.05 and the uniform model with F = 50 Jm“^K“^s“^/^, even 
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Temperature vs. Solar elevation angle 



Solar elevation angle (deg) 



Solor elevation ongle (deg) 


Fig. 15. Same as Fig. 10 but for the core model. In the upper and lower 
panels, L/R = 0.2 and 0.05, respectively. The results of three different radii 
R = 0.3(solid), l.O(dotted), and 3.0 m (dashed) are shown for the southern 
(blue) and northern (red) faces. 
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Fig. 16. Best-fit parameters for the core model. Top: particle radius R. 
Middle: ratio of regolith thickness to particle radius, L/R. Bottom: goodness 
of hts, XBrnin/ximim relative to that of the uniform model. In the middle 
panel, the mean internal density is also shown on the right axis. The ranges 
of R and L/R are limited a priori in data hts as 0.125 m < i? < 8 m and 
0.05 < L/R < 0.5. 
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Fig. 17. Data fitting with the core model at 125,000 km. Top: best-ht 
temperature curves as a function of B'. The northern and southern face 
temperatures are represented by red and blue colors, respectively. Bottom: 
contour of Xb/Xb min~ ^ plane of L/R vs. R. The point of the best-ht 

L/R and R is shown by the plus mark. 
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Fig. 18. Same as Fig. 17, but for the case of a = 129, 000 km 
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though the modeled equinox temperatures are similar to each other. In the 
core model, we see a sharp kink at the equinox, particularly for R = 3 m, 
as the low F regolith responds quickly to the incoming flux. On the other 
hand, the temperature minimum comes slightly after the equinox and the 
temperature curve there is round, if F is large in the uniform model. 

The best-£t parameters for the core model are shown in Fig. 16. The 
particle radius is about 1 m across throughout the A ring, and is slightly 
smaller than those estimated in the uniform model on average. Unlike the 
uniform model, the core model gives upper limits of R in the middle A ring, 
where relatively large values of F are seen in the uniform model. The regolith 
thickness L relative to R takes small values in the middle A ring whereas it 
is high in the inner and outer A ring. In principle, a small L/R is seen at 
the location with a high F in the uniform model. The mean p of the particle 
calculated from L/R is shown on the right axis of the middle panel. The 
mean p is a maximum ~ 800 kg m“^ in the middle A ring and much lower 
values of p (< 300 kg m“^) are seen near the edges of the A ring. The value 
of the mean p in the middle A ring is probably too high to prevent ring 
particles from gravitational accretion, and a finite porosity seems necessary 
for the cores. The bottom panel of Fig. 16 shows XB.min relative to ya min . 
This is about unity at all radial locations, indicating the model fits are good 
also for the core model. 

Examples of best-fit curves and the contours of min — 1 on the L/R 

vs. R plane are shown in Fig. 17 (a = 125, 000 km) and Fig. 18 (a = 129, 000 
km). At a = 125, 000 km, the best-fit parameters are L/R = 0.16 and i? = 91 
cm. Although the error bars of individual parameters are large, we can ex¬ 
clude a combination of a large R and a low L/R] this case gives an equinox 
temperature that is too large compared with the observed temperature. The 
lower limit of R is well constrained as a small R gives a low equinox tem¬ 
perature. At a = 129,000 km, the best-fit parameters are L/R = 0.06 and 
R = 109 cm. The contour of Fig. 18 is similar to that in Fig. 17, but a band 
of low is relatively shifted toward the upper-left corner. The upper limits 
of R and L/R are well constrained in this case, as shown in Fig. 16. It is 
fonnd that if R is large, the core continnously warms up the regolith mantle 
even long after the eqninox. As a result, regardless of L/R, the modeled tem¬ 
peratures, particularly for the south face, become higher than the observed 
temperatnres at B' ~ 2.5° (see the curves of i? = 3 m in Fig. 15). If R is 
small (< Im) this effect is weak and the lower post-eqninox temperatnres are 
obtained. 
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Fig. 19. Median particle sizes deduced from stellar and radio occupations 
(colored diamonds). The horizontal error bar is the range of saturnocentric 
distance used for averaging of observational data or fitted parameters. For 
Cuzzi et al. (2009), we use Rmin = 3 mm, i?max = 10 m, and q = 2.7 for 
the inner A ring and i?min = 4 mm, i?max = 3 m, and q = 3.2 for the outer 
A ring (their Fig. 15.3). For French and Nicholson (2000), q = 2.9 is used 
for the outer A ring, while the error bar for i?med comes from a possible 
uncertainty in q (2.75 < q < 3.1), as discussed by them. For Harbison et 
al. (2013), q = 2.8 is used, while the the error bar for R^ed comes from a 
possible uncertainty in q (we adopt 2.7 < q < 2.9); i?min is altered with q 
using their Fig. 16. The black asterisks are the particle sizes estimated from 
our core model (Fig. 16). 

5. Discussion 

5.1. Ring particle size 

The best-fit particle size is about 1 m in our models, although relatively 
smaller sizes are possible near the inner and outer edges of the A ring. The 
particle size has been estimated from various occultations (Zebker et al., 
1985; French and Nicholson, 2000; Cuzzi et al., 2009; Harbison et al. 2013). 
These studies showed that ring particles have a size distribution with an 
differential power-law exponent of ~ —3, while our model adopts a single 
representative size. Since the thermal flux is proportional to the total cross 
section of particles in the absence of mutual shading, it may be meaningful to 
use the median radius of particles i?med for comparison, where i?med is defined 
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so that the total cross section of particles with R > Rmed is equal to that of 
R < -Rmed- The median radius is given as where 

/?max and -Rmin are the upper and lower cut off radii, and —q is the power- 
law exponent for the differential size distribution. The values of i?med from 
occultations are shown in Fig. 19 together with the particle sizes estimated 
by our core model (Fig. 16). 

For the region inside the Encke division, i?med from all studies are com¬ 
patible with our estimates. For the outer A ring, i?med from French and 
Nicholson (2000) is compatible with our estimates whereas i?med from Cuzzi 
et ah (2009) and Zebker et al. (1985) is much smaller than the lower limit of 
R estimated by our models. The difference in i?med between these studies of 
occultations comes from the population of particles smaller than 1 cm; i?min 
= 1 cm in French and Nicholson (2000) while i?min = 4 mm and 1mm in 
Cuzzi et al. (2009) and Zebker et al. (1985), respectively. 

If the cross section of the outer A ring is dominated by mm to cm-sized 
particles, as suggested by Cuzzi et al. (2009) and Zebker et al. (1985), the 
ring should cool down close to the equilibrium temperature at the equinox. 
On the other hand, the CIRS data show the ring temperatures slightly higher 
than the equilibrium equinox temperatures in the outer A ring (by 3K). 
This may indicate that our radiative transfer models underestimate the equi¬ 
librium equinox temperatures, which we do not believe very likely (Section 3). 
The apparent contradiction may be explained if we consider possible gravita¬ 
tional aggregates (e.g., Salo, 1995) whose total cross section is small but total 
mass dominates in the outer A ring. They can be much warmer than the 
equilibrium temperature at the equinox. If they break up due to mutual col¬ 
lisions (Karjalainen, 2007) for example, small warm particles are dispersed. 
If some of aggregates are persistent on a time scale longer than the thermal 
relaxation time (~ a year), the warm equinox temperature may be explained. 

5.2. Radial variation of internal density 

The temperature anomaly at the equinox is prominent in the middle A 
ring and this was interpreted by a relatively large thermal inertia for the 
uniform model and a relatively large core size in the core model. As pointed 
out in Section 4.5.1, the uniform model is not self-consistent as a very small 
K estimated by the model is not expected for a medium with a moderately 
high density suggested from self-gravity wakes. Therefore, we prefer the 
interpretation of the core model to that of the uniform model. Namely, the 
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internal density of ring particles, p, is highest in the middle A ring and it 
decreases towards the inner and outer edges of the A ring. 

Some other observations seem to support this idea. The amplitude of 
azimuthal brightness asymmetry in photometric observations has a peak in 
the middle A ring (Dones et ah, 1993; French et ah, 2007) and it sharply de¬ 
creases towards the inner and outer edges of the A ring. French et al. (2007) 
suggested several possibilities to produce the relatively weak asymmetric fea¬ 
ture in the inner A ring as compared with the middle A ring: p is lower, the 
particle size distribution is broader, or the restitution coefficient is higher in 
the inner A ring. Among them, any clear difference in the size distributions 
between the inner and middle regions of the A ring is not identihed (Zebker et 
ah, 1985; French and Nicholson, 2000). Also, similar elastic properties of ring 
particles are expected in these regions because similar surface regolith prop¬ 
erties are suggested from ring spectra (e.g., Hedman et ah, 2013). Therefore, 
the radial p-variation seems to be the most plausible cause of the steep radial 
change of the azimuthal brightness asymmetry. Low-p particles in the inner 
A ring are also indicated from existence of axisymmetric wakes (Thomson et 
ah, 2007; Hedman et al. 2014), which are probably produced by overstabil¬ 
ity that is induced if p is low (Salo et ah, 2001). It is unclear whether low-p 
particles are necessary for the weak azimuthal brightness asymmetry seen 
in the outer A ring because other explanations may work, such as a broad 
size distribution, particle clumping rather than wake formation, and velocity 
enhancement due to satellite perturbations. 

Cassini ISS observations identihed many propeller-forming moonlets in 
the A ring (Tiscareno et ah, 2006, 2008, 2010; Sremcevic et ah, 2007). They 
are particularly clustered in the middle A ring (Sremcevic et ah, 2007; Tis¬ 
careno et ah, 2008) that is around the equinox temperature peak. These 
moonlets probably have dense cores surrounded by relatively less dense par¬ 
ticles as suggested for small inner satellites (Porco et ah, 2007; Yasui et ah, 
2014). It is natural that dense particles are also relatively abundant in the 
background of the middle A ring. 

The question is what causes the radial variation of p or regolith thickness, 
if it is real. Unless some mechanisms work, the radial variation of p must 
be smoothed out on the viscous diffusion time, ~ 10® years for the A ring 
(Esposito, 1983; Tiscareno et ah, 2007). Some explanations can be consid¬ 
ered for the difference between the middle and outer regions of the A ring. 
Relatively abundant small particles in the trans-Encke region (French and 
Nicholson, 2000, Cuzzi et ah, 2009) indicate that mutual impact velocities 
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are large there probably due to strong self-gravity wakes or satellite pertur¬ 
bations. It is likely that the dust production rate in such an environment is 
high and ring particles are covered by thick regolith. Accumulation of E ring 
particles may also be important near the outer edge of the A ring (Horanyi 
et ah, 2009). 

The more puzzling mystery is the difference between the inner and middle 
regions of the A ring. We speculate that moonlets selectively transport high-p 
particles outward. Moonlets are known to radially migrate due to collisional 
and gravitational interactions with surrounding particles, although the exact 
migration mechanism is under debate (Crida et al. 2010; Pan and Chiang, 
2010; Rein and Papaloizou, 2010; Tiscareno, 2012; Pan et al. 2012; Bromley 
and Kenyon, 2013). We do not discuss the detailed mechanisms of migration, 
but simply consider a situation that moonlets can migrate a long distance 
either inward or outward on a long timescale due to gravitational torques. 
If moonlets are rubble-pile aggregates, however, inward migration may be 
inefficient as they break up due to the tidal force. Such a restriction does not 
exist for outward migration, although moonlets may not be able to migrate 
beyond some strong satellite resonances (Bromley and Kenyon, 2013) and 
are trapped in certain regions which may be currently seen as the moonlet 
belts in the middle A ring. After moonlets migrate outward, some of them 
may break up due to mutual collisions or meteoroid impacts and disperse 
high-p particles. 

An alternative interpretation for the radial p variation is that the A ring, 
or at least the middle part of it, is simply younger than the viscous diffu¬ 
sion time. For example, a large icy satellite (~ 100 km in size) might have 
existed at the present location of the middle A ring and been destroyed by 
a meteoroid impact. The dispersed chunks of dense material have not had 
enough time to diffuse throughout the A ring. A similar discussion is found 
for radial diffusion of moonlets (Crida et ah, 2010). A young A ring is also 
indicated from other various estimates (Esposito, 1986; Charnoz et ah, 2009 
and references therein). 

In addition to them, thin regolith mantles found in the present work 
provide a new constraint on the A ring age. Elliott and Esposito (2011) 
show that the regolith depth of ring particles increases with time due to 
meteoritic bombardment if ring particles are originally regolith-free. If we 
apply the regolith thickness of 10 cm for 1 m particles (Fig. 16) to the time 
evolution of the regolith thickness (their Fig. 2), the ring age is found to 
be only ~ 10^ yr. Since the meteoroid bombardment flux is recently found 
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to be much smaller (Kempf et ah, 2013) than the previous estimate (Cuzzi 
and Estrada, 1998), the timescale may be larger by one to two orders of 
magnitude. This is still much shorter than the age of the solar system and 
is close to that suggested from the viscous diffusion time. 

A problem of the young A ring hypothesis is that the size of an impactor 
that can destroy a 100 km body is ~ 10 km and that the flux of such massive 
bodies is very low; the probability of destruction of a 100 km body in 10® yr 
is only ~ 10“® (Charnoz et al., 2009). In this estimation, the critical impact 
energy necessary to disperse a half target mass in free space is used (Benz and 
Asphaug, 1999). Impact simulations including the tidal force of the planet 
(Karjalainen, 2007; Hyodo and Ohtuski, 2014), however, show that impact 
outcomes at the A ring location are much more erosive than those in free 
space, if the target is a gravitational aggregate. Thus, destruction of a large 
satellite may not be an unrealistically low-probability event. 

5.3. Composition of particle cores 

We assumed rigid water ice for the core material in the core model. How¬ 
ever, similar values of R and L/R can be derived even for other materials, as 
long as pC for the core is similar to the assumed value [see Eq. flTOj) ]. In fact, 
the values of pC of chondrites (or rocks) are similar to that for water ice; 
as compared with water ice, their densities are larger by a factor of 3-4 but 
the specihc heats are lower by a factor of 3-4 (Yomogida and Matsui, 1983; 
Consolmagno et ah, 2013). If the cores of ring particles are chondritic and 
L/R 0.9, the mean density results in ~ 3,000 kg m“®. Such high density 
particles cannot exist as ring particles at the A ring location, because ac¬ 
cretion must quickly occur. Therefore, we can exclude cores predominantly 
made of rocks. The cores must be predominantly water ice, probably with 
hnite porosities. 

5.f. Other rings 

As found in Spilker et al. (2013), the equinox temperatures for the B 
and C rings are well reproduced by the models with the Saturn flux only. 
For the B ring, we also check the asymmetry between the pre- and post¬ 
equinox temperatures, as done for the A ring (Figs. 7 and 8). We hnd that 
the pre-equinox temperatures are higher than the post-equinox temperatures 
for the unlit face, but it is opposite for the lit face, unlike the middle A ring. 
These facts indicate that the equinox temperature anomaly for the B ring is 
much smaller than that for the middle A ring, even if it exists. This means 
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that the B ring particles are either smaller or less dense than the middle A 
ring particles. Since the particle sizes for the B ring estimated from stellar 
and radio occultations (French and Nicholson, 2000; Cuzzi et ah, 2009) are 
almost the same with A ring particles inside the Encke division, it is likely 
that the internal density of B ring particles is lower than that for middle 
A ring particles or may be similar to inner A ring particles. It is difficult 
to say anything about the temperature asymmetry for the C ring because 
the equinox temperatures strongly depend on observational geometry. The 
correction of the geometry dependence remains for future work. 

6. Conclusion 

In the present paper, we examined the seasonal temperature variation of 
Saturn’s A ring. The lowest temperature was observed at the solar equinox 
in August 2009. We used two different radiative transfer models (a multi¬ 
layer model and a mono-layer wake model) and found that both models give 
the equinox temperatures that are much lower than the observed equinox 
temperatures, as long as only the flux from Saturn is taken into account. We 
also found that the post-equinox temperatures are lower than the pre-equinox 
temperatures at the same absolute solar elevation angle. These facts strongly 
indicate that the A ring was not completely cooled down at the equinox due 
to the effect of a hnite seasonal thermal inertia. 

Since the seasonal thermal skin depth is comparable to the ring particle 
size, it is possible to give constraints on the physical properties of particles 
down to the deep interior using the seasonal temperature variation curves. 
We developed a simple model for seasonal temperature variation, and time 
dependence was taken into account by solving the thermal diffusion equation 
for the temperature evolution of the particle interior. The model calculates 
only the lit and unlit face temperatures, ignoring geometry dependence and 
diurnal variation. We considered two types of the internal structure of ring 
particles. In the first model, uniform internal density and thermal inertia 
of a ring particle were assumed. The particle size was estimated to be 1-2 
m. The seasonal thermal inertia was found to be 30-50 Jm ^s in 
the middle A ring whereas it is ~ 10 Jm“^K“^s“^/^ or as low as the diurnal 
thermal inertia in the inner and outermost regions of the A ring. In the 
second model, a ring particle has a high density core surrounded by a much 
less dense regolith mantle. This model showed that the particle size is still 
about 1 m and that the core radius relative to the particle radius is about 
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0.9 for the middle A ring and is much less for the inner and outer A ring. 
This means that a radial variation of the internal density of ring particles 
exists across the A ring. 

To conhne high density particles in the middle A ring against viscous 
diffusion, some mechanisms may be working. We speculated that moonlets 
are transporting dense particles to the limited regions. An alternative idea is 
that the A ring was recently formed (< 10® yr) by catastrophic destruction 
of a pre-existing icy satellite, so that dense particles have not yet diffused 
over the A ring and regolith mantles of particles have not grown thick. We 
also excluded cores predominantly made of rocks, as such particles have very 
large internal densities. 

Acknowlegements 

We thank anonymous reviewers for helpful comments that signihcantly 
improved the manuscript. This research was carried out at the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with NASA. 
Government sponsorship acknowledged. This work was also supported in 
part by the NASA’s OPR Program (NNX14A036G). The numerical simula¬ 
tions were performed using the JPL supercomputer, Aurora. 

Appendix A 

In the seasonal model introduced in Section 4, the heliocentric distance 
rs and the solar elevation angle B' as a function of time, t, are necessary. We 
simply assume that Saturn has a Keplerian orbit around the Sun and adopt 
the orbital elements at the epoch J2OOO0; the semimajor axis of Saturn as is 
9.537 AU and the orbital eccentricity es is 0.054. The heliocentric distance 
rs is given as (Murray and Dermott, 1999) 


rs = as(l - es cosEs), 


(16) 


where Es is the eccentric anomaly of Saturn. The eccentric anomaly Es as 
a function of time is given by solving the Kepler’s equation: 


^^s(^ - tph) = Es-es sin Es, 


(17) 


Tttp: / /nssdc. gsfc. nasa.gov/planetary / factsheet/saturnfact. html 
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where us = 27r/Ts is the mean motion of Saturn (where Ts = 10759.22 day 
is the orbital period of Saturn) and tph is the time at the perihelion passage. 
Since rs =9.433 AU and Saturn was receding from the Sun at the equinox, 
Es at the equinox was 1.367. We dehne ust = tt at the equinox, and that 
gives custph = 1.827. The angle between the vector pointing towards Saturn 
from the Sun and the Saturn spin vector is dehned to be 7 . Using 7 , the 
solar elevation angle B' is dehned as 5' = 7 — 7 r /2 and given as 

sin B' = sin B'^^^ sin (/g - fs,i) , (18) 


where B'^^^ = 26.73° is the obliquity of Saturn or the maximum solar eleva¬ 
tion angle. In the above, /g is the true anomaly of Saturn given as 


cos/g 


cos Es — eg 
1 — eg cos Es ’ 


(19) 


and /g_i = 1.421 is the value at the equinox. Note that B' increases during 
the half orbit of Saturn in which Cassini observations have been made. 
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